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Abstract 

The form of the initial value constraints in Ashtekar's hamiltonian 
formulation of general relativity is recalled, and the problem of solv- 
■ ing them is compared with that in the traditional metric variables. 

CN . It is shown how the general solution of the four diffeomorphism con- 

straints can be obtained algebraically provided the curvature is non- 
' degenerate, and the form of the remaining (Gauss law) constraints is 

discussed. The method is extended to cover the case when matter 
is included, using an approach due to Thiemann. The application of 
the method to vacuum Bianchi models is given. The paper concludes 
with a brief discussion of alternative approaches to the initial value 
(3JT)! problem in the Ashtekar formulation. 



1 Introduction 



It is with great pleasure that we dedicate this paper to Dieter Brill, our 
teacher, advisor, and colleague, on the occasion of his 60th birthday. Our 
contribution concerns the initial value problem for general relativity, which 
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is amongst Dieter's many areas of expertise. As is the case with almost all 
research activity developed around the general relativity group at the Uni- 
versity of Maryland, the ideas we will present have benefitted from Dieter's 
always kind and sometimes maddening insightful questioning. Of course it 
is our wish that this paper will prompt some more such questioning. 

General relativity is invariant under four dimensional diffeomorphisms. In 
the hamiltonian formulation, this invariance manifests itself in the presence 
of constraints on the canonical variables. If the constraints are satisfied at 
one instant of time, they continue to hold at all times. The initial value 
problem is the problem of giving a construction for a parametrization of the 
general solution of the constraints. 

The most well-developed approach to the initial value problem is based on 
the "conformal technique" (see for example |1J and references therein) . In this 
approach, the phase space variables for general relativity are the 3-metric q a b 
and its conjugate momentum p ab which is a tensor density of weight 1, closely 
related to the extrinsic curvature of the spatial hypersurface in a solution to 
the field equations. There are four constraints, consisting of one Hamiltonian 
constraint, which generates diffeomorphisms normal to the spacelike surface 
on which the initial data are given, and three momentum constraints, which 
generate spatial diffeomorphisms. The Hamiltonian constraint is viewed as 
a (quasi-linear, elliptic) equation determining the conformal factor of the 
3-metric. The momentum constraints determine the "longitudinal" piece of 
p ab . In this approach, the freely specified data are, in principle, the conformal 
equivalence class of the 3-metric, and the transverse traceless part of the 
momentum p ab . 

The conformal technique is well suited to addressing questions of existence 
and uniqueness of solutions to the field equations. However, for two reasons 
it is not suited to reduction of the phase space to the unconstrained degrees 
of freedom. First, the "freely specifiable" transverse traceless tensor densities 
are not known for an arbitrary conformal 3-metric. Second, the solutions to 
the constraint equations determining the conformal factor and longitudinal 
part of the momentum are not given explicitly. 

An alternative representation for hamiltonian general relativity was intro- 
duced in the mid-eighties by Ashtekar ||. The Ashtekar representation may 
be understood (and was discovered originally) as resulting from a complex 
canonical transformation that goes from the (triad-extended) ADM gravi- 
tational phase space variables to a new set of variables. The new variables 
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are a spatial S0(3,C) connection, and an so(3,C)-valued vector density as 
its conjugate momentum. The major benefit of this representation is that 
the constraints take a simpler form in terms of these variables than in the 
ADM representation. First, the constraints turn out naturally to be polyno- 
mial of low order in the Ashtekar phase space variables. Moreover, for the 
vacuum case, the scalar constraint is homogenous in the Ashtekar canoni- 
cal momenta. This simplification has allowed considerable progress in many 
questions of gravitational physics. Perhaps the most impressive results have 
been obtained in the quantum theory, where this formulation has made it 
possible to perform the first few steps in the Dirac quantization of general 
relativity §@,g. 

Given that the Ashtekar representation simplifies the form of the con- 
straints, it is natural to reconsider the classical initial value problem in this 
context. As it turns out, generically in phase space, it is possible to obtain 
algebraically the general solution to the Ashtekar version of the diffeomor- 
phism constraints of general relativity! || |7j However, this does not give an 
easy solution of the initial value problem. The reason is that, due to the 
covariance of the Ashtekar formalism under local SO(3,C) rotations, there 
are three additional constraints whose form is identical to the non-abelian 
form of Gauss' law familiar from Yang-Mills theories. Once the diffeomor- 
phism contraints are solved by the method to be described below, one still 
has to face the Gauss law, which has now become a non-linear first order 
partial differential equation on the remaining variables. In addition, the re- 
ality conditions restricting to the phase space of real, Lorentzian GR need to 
be enforced. 

Still, one can hope that the relocation of the difficulties afforded by the 
Ashtekar approach may prove to be fruitful in some applications. This al- 
ready seems to be the case in the context of GR reduced by symmetry con- 
ditions. In addition, the interplay between "gauge-fixing" and solving the 
constraints is different in the Ashtekar approach, which is a fact that re- 
mains to be fully explored. 

The rest of this paper is organized as follows. The form of the initial value 
constraints in Ashtekar's hamiltonian formulation of general relativity is re- 
called in section 2, and in section 3 it is shown how the general solution of the 
four diffeomorphism constraints in vacuum (or with a cosmological constant) 
can be obtained algebraically provided the curvature is non-degenerate. The 
form of the remaining (Gauss law) constraints is also discussed, as is the 
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relation between the solution given and the structure of the 4- dimensional 
curvature. In section 4 the method is extended to cover the case when mat- 
ter is included, using an approach due to Thiemann. The application of the 
method to vacuum Bianchi models is given in section 5, and the paper con- 
cludes in section 6 with a brief discussion of alternative approaches to the 
initial value problem in the Ashtekar formulation. 

2 Ashtekar 's variables 

Ashtekar's representation of hamiltonian general relativity has been the sub- 
ject of extensive reviews (cf. |, |), so we shall just recall its main features. 

The Ashtekar canonical coordinates are an SO(3,C) spatial connection, 
A l a , and an so(3,C)-valued vector density of weight 1, Ef. The fundamental 
Poisson bracket is given by {A l a (x), Ej(y)} = iSj6^P(x,y). (Our notation is 
as follows. Latin letters from the beginning of the alphabet denote spatial 
indices, e.g. a,b,... = 1,2,3. Latin letters from the middle of the alphabet 
are SO(3,C) indices, = 1,2,3. They are raised and lowered with the 

Kronecker delta and 6ij. We will also use the totally antisymmetric symbol 
6ijk, with ei23 = 1. We use units with c = G — 1.) 

The vector density E® may be identified with the densitized spatial triad 
y^e*, with the contravariant spatial metric given by q ab = e"e 6 \ In turn, the 
SO(3,C) connection may be identified in a solution with the spatial pullback 
of the self-dual part of the spin-connection. 

These variables parametrize the phase space of complex general relativ- 
ity. A real metric with Lorentzian signature may be recovered by imposing 
appropriate reality conditions (cf. ||). These conditions amount to the re- 
quirement that E^E hl be real, and that its time derivative be real. If these 
conditions are satisfied initially then the dynamical evolution will preserve 
them in time. A metric of Euclidean signature is obtained by simply taking 
A\ and Ef real. 

In terms of these variables, the constraints of (complex) general relativity 
take the form 



l3k E a t E)B c k 

r- T?b TDci 
DaE? 



—iJ, 



a 



(1) 

(2) 
(3) 
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Here e a bc is the standard Levi-Civita tensor density. Bf is the "magnetic 
field" of the connection A l a , defined by Bf := e abc F bci , and F l ab := d a A\ — 
dbA l a + e l jkA J a A b . D a is the SO(3)-covariant derivative determined by A % a . 
p/y/Q an d Ja are (up to coefficients) the matter energy and momentum den- 
sities respectively. Ki is a spin density, present only for half-integer spin 
matter fields. 

In the following, we will call the constraints ([!]) and (§) "diffeomorphism 
constraints" . The first generates diffeomorphisms normal to the spatial hy- 
persurface S together with some SO (3) rotation, so it takes the place in 
Ashtekar's formalism of the ADM hamiltonian constraint. It will be called 
here the "scalar constraint". The second generates diffeomorphisms tangen- 
tial to E together with some SO (3) rotation, so it takes the place of the ADM 
momentum constraint. It will be called here the "vector constraint". 

There are 3 additional constraints, ([|), of the Gauss- law type familiar 
from Yang-Mills theories. They generate local SO(3,C) rotations, under 
which the Ashtekar formalism is covariant. 

The constraints are polynomial in the gravitational phase space variables 
and for scalar and spin- 1/2 matter. For Yang-Mills type fields, polynomiality 
is mantained only if one multiplies the scalar constraint ([!]) through by detq = 
detEf, which is cubic in Ef. Moreover, it is remarkable that the gravitational 
contribution to the scalar constraint (|1]) is homogenous in the momenta Ef. 
This should be compared with the ADM Hamiltonian constraint, where the 
term quadratic in the canonical momenta must be balanced by a "potential" 
term given by the 3D Ricci scalar times the determinant of the 3-metric. 
Note that while the vector and Gauss constraints @) and fl3|) are densities of 
weight 1, the scalar constraint (|l|) is of weight 2. 

3 Algebraic solution of the vacuum diffeo- 
morphism constraints 

We now proceed to show how one can use the simplification of the constraints 
for general relativity provided by the Ashtekar formalism to tackle the ini- 
tial value problem. In particular, we will show how the general solution of 
the the scalar and vector diffeomorphism constraints can be obtained by al- 
gebraic methods. We originally discovered this solution in the context of 
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a Lagrangian pure spin-connection formulation of GR, in which one solves 
for the metric variables (self-dual 2- forms) in terms of the connection |J. 
However, it is unnecessary to view the technique in that context. 

The first step is to assume that the magnetic field Bf is non-degenerate 
as a 3 x 3 matrix. In a generic real, Lorentzian spacetime, the real and 
imaginary parts of the complex equation detB = will define two surfaces, 
and their intersection will give a one-dimensional submanifold on which Bf 
is degenerate. That is, generically Bf is non-degenerate except on a set 
of measure zero. (The points in phase space with detB = everywhere 
form a set of measure zero in phase space. The identity of these points is 
best seen in the covariant formalism (cf. || It turns out that they 
correspond to space-time metrics of Petrov type 0, (4), (3,1), when Ef is 
non-degenerate.) We will ignore here any problems that might arise as a 
consequence of degeneracy of Bf. 

Now, since Bf is assumed to be non degenerate, one may use it as a basis 
for the space of SO(3,C) vector densities. In particular, one can expand the 
momentum Ef with respect to it, 

Ef = P ld B a > , (4) 

for some (non-degenerate) 3x3 matrix Pij. 

We first consider the vacuum case, p = 0, J a = 0, Ki = 0. It is easy to see 
that the vector constraint @ implies that must be symmetric. The scalar 
constraint ([!]) implies that P^ must satisfy a quadratic algebraic condition, 

(Pf) 2 - PijP 13 = . (5) 

This condition fixes one of the six components of P^ with respect to the 
others. 

We know of two methods for solving the scalar constraint (|5|). In the first 
method, one solves @ for the trace of P. Decomposing P^ into the trace 
TrP and tracefree part P^ := P^ — \TrP5ij, (|j) becomes the condition 

TrP = ±(|TrP 2 ) 1/2 . (6) 

In the second method for solving ([5]), one notes that the characteristic 
equation for P implies that TrP' 1 = [TrP 2 - (TrP) 2 ]/2detP. Thus (|) 
is equivalent to the tracelessness of P" 1 , provided P is invertible (which 
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it must be if both Bf and the 3- metric are nondegenerate) . Let if> denote 
P -1 . Then the characteristic equation for if), assuming Trip = 0, reads 
ip 3 - \{Trif) 2 )if) - detipl = 0, giving if)~ l = [ip 2 - \(Trif) 2 )I]/ detif). Thus the 
general solution to @ for invertible P can be explicitly written in terms of 
the 5 independent components of a tracefree 3x3 matrix <f> in the form 

P = 2 -i(Tr0 2 )/. (7) 

One is left with the Gauss constraint @, which becomes 

BfDaF** = , (8) 

where we have used the Bianchi identity D a B% = 0. In view of (^) or 
the Gauss constraint becomes a non-linear equation in P^- or respectively. 
It is here, and in the reality conditions, that the initial value problem now 
resides. 

We can consider the Gauss law constraint as 3 conditions on the 5 in- 
dependent components of P % \ for a fixed A % a . Then there are presumably 2 
free functions worth of solutions for this equation, yielding the two indepe- 
nent metric degrees of freedom of GR. Up to this stage the connection A l a 
has remained entirely arbitary. By use of the 4 parameter diffeomorphism 
and 3 parameter SO (3) transformations, one could now fix all but 2 of the 9 
components of the connection, yielding the other half of the coordinates on 
the reduced phase space. 

If a spherically symmetric ansatz is assumed, then the above method of 
solving the constraints can be carried out entirely, and the Gauss constraint 



reduces to an ordinary differential equation that can be solved ||10|| . 

The existence in general of solutions to the Gauss law constraint (§) with 
P* J restricted by @ or (0) is a problem that has not been addressed. To 
indicate just one complication that can arise, consider the asymptotically 
flat case. If A\ is asymptotically the spin-connection corresponding to a 
negative-mass spacetime, then we know by the positive energy theorem that 
it must be impossible to find a regular solution to (^) for P iJ , subject to the 
reality conditions. 

The geometrical interpretation of the matrix P y is easily seen from the 



covariant point of view [II]. The densitized triad E a% of the Ashtekar formal- 
ism is related to the triad of anti-self-dual 2- forms by E ai = e abc |L2 1 . Now 



the curvature 2-form R % of the spin-connection can be expanded in terms of 
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the self-dual and anti-self-dual 2- forms as as Ri = ipij E J + \A Ej + Qfj^ ■ 
is symmetric and tracefree and is just the Weyl spinor in SO (3) notation. A 
is proportional to the Ricci scalar, and is equivalent to the tracefree part 
of the Ricci tensor. Thus in vacuum, A and $jj vanish, and one can solve for 
the ant i- self- dual 2-forms in terms of the curvature as S* = (ifj -1 )^ Rj. The 
dual of the spatial pullback of this equation yields immediately the general 
solution given above for the four diffeomorphism constraints in Ashtekar's 
formalism, with P 13 identified with (ip ) %3 . 

Let us now consider the addition of a cosmological constant A. This 
modifies only the scalar constraint, which becomes 

e ahc ^ k EtE)Bl - {\l?>)ke ahc ^ k EtE)El = . (9) 

One may now proceed as in the vacuum case. The only difference is that the 
algebraic condition (Q) is replaced by 

(Pi) 2 - / , ,,/ , '•' = 2A{detP) . (10) 

This can be regarded as a cubic equation for TrP. It is equivalent to the 
statement that the inverse of P has trace equal to A. In this case, P corre- 
sponds to the inverse of the matrix + ^ASij). 

An interesting special solution of the constraints in the case of a nonva- 
nishing A is given by the Ansatz = (3/A)5ij, or 

E? = (3/A)B? , (11) 

which solves automatically all of the constraints. (The Gauss constraint is 
satisfied as a consequence of the Bianchi identity D a B^ = 0.) This Ansatz 
was first introduced by Ashtekar and Renteln who observed that it gives 



rise to self-dual solutions of the Einstein equation with a non- vanishing A. 



4 Matter couplings 

It is possible to extend the method just described to solve the diffeomorphism 
constraints in the presence of matter using an approach that has recently 
been brought to our attention by Thomas Thiemann 0. Thiemann's method 
proceeds as follows. 
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One begins with the constraints in the form ([!]), (|2]), fl3]) ||14|| . Then ex- 
panding E" 1 = PijB a i as in @, one finds that the vector constraint @ implies 
that the anti-symmetric part of Py no longer vanishes, but is given instead 
by 

Mj '■= P[ij] = 2£ e ijkBl J a , (12) 

where B := detB^. This is an explicit expression for Pm provided J a does 
not depend upon P^-j as well. Since J a is just the generator of spatial diffeo- 
morphisms for the matter variables, it does not depend on the gravitational 
field variables for the case of integer spin matter. In the spin-1/2 case, the 
diffeomorphism is accompanied by an SO (3) rotation, which involves the 
spin connection but not Ef. In all cases therefore, J a is independent of P^ . 
(In particular, in the scalar, spin-1/2, and Yang-Mills cases, the structure 
of J a is ird a (j), ii A D a ipA, and e bI f a bi respectively, where e bI and f a bi are the 
Yang-Mills-electric and (dual of) magnetic fields.) 

Now turning to the scalar constraint, we begin by noting that when P^ 
is decomposed into its symmetric and antisymmetric parts Py = SV,- + Aij, 
the left hand side of the scalar constraint ([!]) is given by 

e ijk e abc E ai E b] B ck = detB [{TrS) 2 - TrS 2 - TrA 2 ] . (13) 

Let us first consider the case of a single scalar field. Then we have 

p = 7T 2 + E ai E\d a ^ + e ljk e abc E ai E b ^E ck V{<j>) . (14) 

Now as in the vacuum case there are two approaches to solving the scalar 
constraint. In the first method, one regards the constraint as a cubic equation 
on the trace of S^, which can be solved (albeit messily) in closed form. 
In the second approach, due to Thiemann 0, one notes that every term 
of the constraint is either independent of the scalar field momentum 7r or 
depends on it quadratically! To see why, observe that J a is linear in it, and 
therefore according to (|i~2"|), so too is A^. Thus the gravitational part of the 
constraint (|T3| ) is quadratic in ir. As for p (|14D, the second (gradient) term is 
independent of n because, as one easily sees using (0), A^ drops out of the 
expression E az d a <f) = S^Bjd a (j). Moreover, in the third term one can show 
that detE involves A^j, and therefore tt, also only quadratically. 
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Thus one can simply solve the constraint for ir in closed form, 
the case V(4>) = 0, the resulting expression is fairly complicated: 



Even for 



7T = ± 



((TrSr-TrS' + (S0 2 )/(^-l) 



1/2 



(15) 



where £* := B ai d a (p. 

For the real theory, there is also the constraint that the argument of 
the square root be a positive real number. When V(<f)) ^ 0, the solution 
for 7r becomes significantly more complicated. It is unlikely that this is 
of any practical use in general. However, Thiemann's intended application 
is quantization of the spherically symmetric scalar-gravitational system, for 
which the resulting expressions may be more useful. 

Note that even when the argument of the square root in (f[5"l) is real, the 
solution for it has a sign ambiguity. This can lead to problems if one wishes 
to pass to the reduced phase space for the purpose of quantization, since it is 
difficult to eliminate the momentum and not lose part of the reduced phase 
space [|3~5j | . 

For other types of matter coupling the situation becomes yet more compli- 
cated. For a massless spin-1/2 field, the scalar constraint remains quadratic 
(but not homogeneous) in the two components n A of the spinor field momen- 
tum, so one can in principle solve for one of these components. In addition, 
the Gauss constraint is augmented by a spinor field term in this case. When 
Yang-Mills fields are included, the scalar constraint remains polynomial only 
if it is multiplied through by detE, resulting in a 4th order polynomial in the 
scalar field momentum. 



5 Bianchi models 

In this section, we apply to spatially homogenous vacuum space-times, i.e. Bianchi 
models, the method for the solution of the constraints illustrated in section 3 
for the vacuum case. The hope is to gain some insight on how to deal with the 
remaining, Gauss constraint in the form (§). In any case, this analysis may 
be of use in a minisuperspace approximation of vacuum general relativity. 

For Bianchi models, as follows from the assumption of spatial homogene- 
ity, this last constraint reduces to an algebraic condition. There are some 
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simplifications, but the condition is still rather complicated. We record it 
here for the benefit of inventory. One interesting feature is that only the 
traceless part of enters in it. 

The formulation of Bianchi models in terms of Ashtekar variables was 



first worked out by Kodama [16|. Here we follow the approach of Ashtekar 



and Pullin 17 



As is familiar from the ADM treatment of Bianchi models (see e.g 
|19||), we consider a kinematical triad of vectors, Xf, which commute with 
the three Killing vectors on the spatial hypersurface S. (Capital latin letters 
I, J,K, ... will be used to label the triad vectors.) The triad satisfies 

[X I ,X J \ a = C IJ K X%, (16) 

where Cu K denote the structure constants of the Bianchi type under con- 
sideration. The basis dual to Xf, xi, satisfies, 

2^ox5 = -Cj^xtx* ■ (17) 

Without loss of generality, one may set 

C I j K = e JJL S LK + 26ffVj ] (18) 

with S IJ symmetric. (This S IJ has nothing to do with the symmetric part 
S l i of P y referred to earlier.) From the Jacobi identities, it follows that 
S IJ Vj = 0. The Bianchi classification may be described in terms of the 
vanishing or not of Vj, and the signature of S IJ , subject to this condition. 
The most popular models are Bianchi I, selected by Cjj k = 0, and Bianchi 
IX, selected by Vi = and S IJ = 5 IJ . 

The Ashtekar gravitational phase space variables may be expanded with 
respect to the kinematical triad Xf, and Xai as 

K = A^x? (19) 
Ef = det X EfX a M (20) 

where detx denotes the determinant of Xai which is introduced in order to 
de-densitize Ef. Similarly, the magnetic field may be expanded as 

B? = det X BfX a M , (21) 
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and Bf is given by 

B? = -e MNP C NP Q A Ql + e MNP e ijk AU k P ■ (22) 

The gravitational phase space has thus been reduced to the matrices {A l M , Ef 1 }, 
i.e. from 18 degrees of freedom per space point to only 18 in total. 
The constraints for vacuum reduce to 

(23) 
(24) 
(25) 

where Bf 1 is given by fl2"2|). 

We can follow now the footsteps of section 3, for the solution of these 
constraints. Assuming that Bf 1 is non degenerate, we expand Ef 1 in terms 
of it: 

Ef = P l3 B M i . (26) 

From (p4|), we find that Pij must be symmetric, and from (p3|) that it must 
satisfy the algebraic condition (|||). We arrive at the last constraint, which 
now takes the form 

C KM K B M 1P 13 + e ijk A j M Bf P lk = . (27) 

Using (j22"D and the Jacobi identities, after some algebraic manipulations, this 
may be written in the form 

W j Pij = Z l]k pi k , (28) 

where corresponds to the traceless part of Pij, and 

W j := 3V Q e QMN e jkl A Mk A Nl (29) 
Zijk '■= —2£ijiA l M S MN A Nk . (30) 

Note that this condition involves only the traceless part of P l K 

We can now go through the Bianchi classification for some understanding 
of this condition. For a Bianchi model of type I, defined with Cu K = 0, this 
condition is trivially satisfied. If S IJ = 0, which defines a model of type V, 
the condition reduces to the requirement that W l be a null eigenvector for 



ijk r?M rpN r>P _ 

tMNP&i -D — 

n K TP M _L ^ A3 T? Mk — 

^KM &i + £ijkA M rLi — 
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Pi4. IfVc 



Q 



0, one can look at the simple case in which the connection matrix 
Am% is assumed to be diagonal (see Ref. jlT]]). Then, recalling that without 
loss of generality S IJ can be put in diagonal form with ±1 or entries, and 
letting S IJ = diag(si, s 2 , s 3 ), the condition takes the form 



[si(A 11 ) 2 -s 2 (A 22 ) 2 }P 12 = 



(31) 



together with its cyclic permutations. An obvious solution is obtained if Pij 
itself is diagonal. 



6 Other Approaches 

The introduction of Ashtekar's variables has generated other approaches to 
the initial value problem in general relativity. One is the recent proposal 



by Newman and Rovelli p0 |. For a Yang-Mills theory they use a Hamilton- 
Jacobi technique to solve Gauss' law. The Yang-Mills physical degrees of 
freedom are coded in one pair of scalar functions per dimension of the gauge 
group, together with conjugate momenta. Each pair of scalar fields defines a 
congruence of lines ( "generalized lines of force" ) by the intersections of their 
level surfaces. The method can also be applied to gravity in the Ashtekar 
formulation, where Newman and Rovelli solve not only the Gauss constraint, 
but also the vector constraint, by the device of using three of the scalar fields 
as spatial coordinates and letting the remaining fields be given as functions 
of these. It is not known at present if also the scalar constraint can be solved 
using this technique. 

Another potential avenue to the initial value problem is the application of 



the Goldstone- Jackiw solution [21] of the Gauss constraint in the SU(2) Yang- 
Mills theory to gravity in the Ashtekar formulation P2[ . In this approach, 
one would start by solving the Gauss law. The solution could then be used 
in the diffeomorphism constraints. For Yang-Mills, this procedure is not 
particularly useful, because it results in a complicated hamiltonian. For 
gravity, the jury is still out. 

Finally, Thiemann |7j has just introduced a new approach that enables 
him to solve all of the constraints. The key idea is to start with the ansatz 
E at = e abc D]pv l c y This ansatz involves no loss of generality provided the 
curvature is non-degenerate, or provided the spin density K { vanishes (so 
that the Gauss constraint implies D a E a% = 0). Using this ansatz, the Gauss 
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constraint becomes a linear condition on v l a . The vector and scalar constraints 
are then solved in the presence of matter by solving for some of the matter 
momenta. 

References 

[1] Choquet-Bruhat Y. and York J.W. Jr., in General Relativity and 
Gravitation, vol. 1, ed. A. Held, Plenum Press, 1980. 

[2] Ashtekar A., Phys. Rev. Lett. 57, 2244 (1986) ; Phys. Rev D 36, 
1587 (1987) 

[3] Jacobson T., and Smolin L. Nucl. Phys. B299, 295 (1988). 

[4] Rovelli C, and Smolin L. Nucl. Phys. B331, 80 (1990). 

[5] Ashtekar A., (in collaboration with R.S. Tate) Lectures on Non- 
Perturbative Canonical Gravity, World ScientiGc, 1991. 

[6] Capovilla R., Dell J., and Jacobson T., Phys. Rev. Lett. 63, 2325 
(1989). 

[7] personal communication, October, 1992; Thiemann T., "On the 
initial value problem for general relativity coupled to matter 
in terms of Ashtekar's constraints", preprint PITHA 93-1, ITP, 
RWTH Aachen, Germany, January 1993. 

[8] Rovelli C, Class. Quant. Grav. 8, 1613 (1991). 

[9] Capovilla R., Dell J., and Jacobson T., Class. Quant. Grav. 8, 59 
(1991); Erratum 9, 1839 (1992). 

[10] Bengtsson L, Class. Quant. Grav. 8, 1847 (1991). 

[11] Capovilla R., Dell J., Jacobson T., and Mason L., Class. Quant. 
Grav. 8, 41 (1991). 

[12] Mason L., and Frauendiener J., in Twistors in Mathematics and 
Physics, eds. Baston R., and Bailey T., Cambridge University 
Press, 1990. 



14 



[13] Ashtekar A. and Renteln P., in Mathematics and General Relativ- 
ity, ed. J. Isenberg, American Mathematical Society, 1988. 

[14] Ashtekar A., Romano J.D., and Tate R.S., Phys. Rev. D 40, 2572 
(1989); Jacobson T., Class. Quant. Grav. 5, L143 (1988). 

[15] The importance of this point was stressed to us by R. Tate in a 
private communication, January 1993. 

[16] Kodama H., Prog. Theor. Phys. 80, 1024 (1988). 

[17] Ashtekar A. and Pullin J., Proc. Roy. Soc. Israel 9, 65 (1990). 

[18] MacCallum M.A.H., in General Relativity. An Einstein Centenary 
Survey, eds. S.W. Hawking and W. Israel, Cambridge University 
Press, 1979. 

[19] Ryan M.P. Jr., and Shepley L.C, Homogeneous Relativistic Cos- 
mologies, Princeton University Press, 1975. 

[20] Newman E.T. and Rovelli C, Phys. Rev. Lett. 69, 1300 (1992). 

[21] Goldstone J. and Jackiw R., Phys. Lett. B 74, 81 (1978). 

[22] Manojlovic N., Syracuse University Preprint, to appear. 



15 



